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We  review  both  the  theoretical  electromagnetic  response  and  the  spectroscopic  measurements  of  metamateri¬ 
als.  To  critically  examine  published  results  for  metamaterial  structures  operating  in  the  range  from  terahertz 
to  optical  frequencies,  we  focus  on  protocols  allowing  one  to  extract  the  optical  constants  from  experimental 
observables.  We  discuss  the  complexity  of  this  task  when  applied  to  metamaterials  exhibiting  electric,  mag¬ 
netic,  and  magneto-optical  response.  The  general  theory  of  the  electromagnetic  response  of  such  systems  is 
presented  and  methods  are  described.  Finally,  we  briefly  overview  possible  solutions  for  implementing 
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1.  INTRODUCTION 

There  has  been  unprecedented  activity  and  excitement  re¬ 
cently  over  the  discovery  of  negative  index  materials 
(NIMs).1  Although  discovered  only  five  years  ago,  re¬ 
search  in  the  field  has  rapidly  blossomed  engaging  a  com¬ 
munity  of  researchers  with  backgrounds  in  electromag¬ 
netics  (EMs),  optics,  engineering,  chemistry,  physics,  and 
materials  science.  Excitement  in  NIM  stems  from  the 
ability  of  these  materials  to  not  only  span  many  decades 
of  the  EM  spectrum,  but  also  to  attain  material  response 
that  is  difficult  or  impossible  to  achieve  with  natural  ma¬ 
terials.  The  additional  inherent  constructed  nature  of 
these  artificial  materials  ensures  applicability  for  future 
devices.2 

The  electrodynamics  of  NIMs  was  first  analyzed  in 
1967  by  Veselago,3  and  the  reversal  of  phenomena  such  as 
refraction,  the  Doppler  effect,  Cerenkov  radiation,  and  ra¬ 
diation  pressure  were  predicted.  Discussion  of  negative 
index  (NI)  was  also  presented  in  the  lecture  notes  of  Man¬ 
delshtam  in  the  early  1940s.4  The  first  experimental  dem¬ 
onstration  of  NI1  was  motivated  by  the  theoretical  work  of 
Pendry  et  al.  who  showed  that  artificially  constructed  ma¬ 
terials  could  be  fabricated  to  exhibit  an  effective  permit¬ 
tivity  and  permeability.5,6 

Metamaterials  constructed  for  NI  behavior  are  typi¬ 
cally  formed  from  individual  components  designed  to  re¬ 
spond  independently  to  the  electric  and  magnetic  fields  of 
an  incident  EM  wave.  The  requirements  for  a  material  to 
obtain  NI  (given  by  the  real  part  of  the  complex  refractive 

index  h  =  nt  4  in.  2)  are  that  the  group  Vg  and  phase  Vp  ve¬ 
locities  be  antiparallel7;  hence  the  index  is  negative,  ri\ 
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<0.  One  such  way  for  a  material  to  meet  this  require¬ 
ment  is  through  negative  values  of  the  real  parts  of  both 
the  electric  permittivity  e1  (of)  and  the  magnetic  perme¬ 
ability  p.i(w)  over  a  common  range  of  frequencies.  Materi¬ 
als  that  exhibit  negative  permittivity  are  not  uncommon, 
and  any  metal  below  its  plasma  frequency  iop  meets  this 
requirement.  Apart  from  metals,  composite  nonconduct¬ 
ing  materials  termed  artificial  dielectrics  have  been 
known  since  the  1940s,  and  are  capable  of  yielding  e<  0 
over  a  limited  frequency  range. 

The  key  to  demonstration  of  NI  was  the  fabrication  of 
an  artificial  material  that  satisfied  the  p.<  0  criterion  in 
simultaneity  with  e<  0.  Although  there  are  natural  mate¬ 
rials  capable  of  negative  magnetic  response  [for  example, 
antiferromagnets  close  to  resonance  frequencies  in  the  far 
infrared  reveal  /x<0  (Ref.  8)],  these  systems  are  usually 
insulating  implying  a  positive  dielectric  function  in  the 
region  where  the  permeability  is  negative.  The  benefit  of 
a  composite  material  is  that  both  response  functions  (e,/r) 
can  be  designed  to  be  negative,  and  thus  one  is  not  con¬ 
strained  by  the  limitations  of  mother  nature. 

The  first  experiments  on  NIM  were  performed  at  micro- 
wave  frequencies  due  to  the  simplicity  of  fabrication  as 
well  as  ease  of  complex  transmission  and  reflection 
experiments.1,9  These  experiments  included  demonstra¬ 
tion  of  one-dimensional  NI,  two-dimensional  NI,  and  veri¬ 
fication  of  negative  refraction.  Since  the  demonstration  of 
NI  there  has  been  a  rapid  extension  of  the  magnetic  com¬ 
ponent  to  both  higher  and  lower  frequencies.  Research 
groups  have  been  pushing  artificial  magnetic  materials  to 
RFs  for  use  in  magnetic  resonance  imaging,10  while 
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others  have  extended  to  terahertz  (THz)11  and  mid- 
infrared  (MIR)12’13  frequencies.  Even  though  the  field  of 
NI  systems  is  still  in  its  infancy,  resonant  behavior  of 
structures  designed  to  exhibit  artificial  magnetic  re¬ 
sponse  have  already  been  shown  to  scale  over  6  orders  of 
magnitude  in  frequency,  thus  evincing  their  usefulness 
for  device  applications.  In  this  paper  we  overview  the 
fundamental  electrodynamics  of  these  resonant  struc¬ 
tures  and  discuss  a  series  of  experiments  aimed  at  the 
demonstration  of  magnetic  response  at  THz  and  higher 
frequencies. 


2.  EXPERIMENTAL  TECHNIQUES  AND 
OBSERVABLES 

At  frequencies  above  the  microwave  regime  experimental 
observables  suitable  for  the  analysis  of  the  EM  response 
of  a  material  include  reflectance  R(co),  transmission  T(oj), 
and  ellipsometric  coefficients  'l'(o))  and  A (co).  The  vast  ma¬ 
jority  of  natural  materials  do  not  show  magnetic  response 
(ji- 1)  and  therefore  the  only  relevant  material  parameter 
is  complex  e.  The  real  and  imaginary  parts  of  e  can  be  in¬ 
ferred  through  one  or  several  of  the  following  procedures: 

(i)  a  combination  ofi?(m)  and  T(co)  obtained  for  trans¬ 
parent  materials  can  be  used  to  extract  the  dielectric 
function  through  analytic  expressions, 

(ii)  Kramers-Kronig  analysis  of  R(to)  for  opaque  sys¬ 
tems  or  of  T(a>)  for  a  transparent  system,14 

(iii)  ellipsometric  coefficients  'I'(oj)  and  A(to)  for  either 
transparent  or  opaque  material  can  be  used  to  determine 
the  dielectric  function  through  analytic  expressions, 

(iv)  various  interferometric  approaches,  in  particular 
Mach-Zehnder  interferometry,15  and 

(v)  THz  time-domain  spectroscopy  (TDS).16 

Apart  from  protocols  (i)-(iii),  which  are  based  on  intensity 
measurements,  the  last  two  techniques  are  capable  of 
yielding  the  optical  constants  directly  through  the  analy¬ 
sis  of  the  phase  information.  While  quite  a  powerful  tech¬ 
nique,  phase  information  for  both  the  reflected  and  the 
transmitted  waves  is  often  difficult  to  obtain  at  higher 
frequencies.  Using  interferometric  methods  it  is  often  pos¬ 
sible  to  obtain  the  phase  of  the  transmitted  wave;  when 
the  sample  is  sufficiently  well  matched  to  free  space,  such 
that  multiple  reflection  within  the  sample  is  minimized, 
then  the  transmitted  phase  is  proportional  to  the  refrac¬ 
tive  index  of  the  sample.  For  matched  samples,  then,  the 
index  can  be  deduced  by  observation  of  the  phase  varia¬ 
tion  with  either  frequency  or  sample  thickness.  This  tech¬ 
nique  has  been  explored  for  NI  media  both  at  microwave 
frequencies17,18  and  in  the  IR.19 

It  is  worth  noting  that  with  the  exception  of  (ii),  all 
these  protocols  produce  the  optical  constants  throughout 
the  entire  frequency  range  where  experimental  data  ex¬ 
ist.  The  standard  routines  (i)-(v)  have  to  be  modified  for  a 
NIM  since  in  this  case  at  least  four  unknown  functions, 
complex  e  and  complex  1,  are  jointly  responsible  for 
the  behavior  of  experimental  observables.  Last  we  point 
out  that  if  one  is  interested  in  characterizing  a  NI  mate¬ 
rial,  a  refraction  experiment  is  suitable  (perhaps  pre¬ 
ferred)  as  this  demonstrates  NI  directly.  Indeed  the  issue 


of  NI  was  largely  debated  even  after  initial  refraction  ex¬ 
periments  were  carried  out,9  until  verified  by  other  vari- 
20 

ous  groups. 

In  Fig.  1  we  show  a  reflectance  spectrum  for  a  Lorentz 
oscillator  modeled  using  the  Fresnel  equations  shown  in 
Eq.  (2),  with  (9=0  and  some  corresponding  optical  con¬ 
stants.  Reflectance  (and  transmission)  is  not  a  particu¬ 
larly  intuitive  material  parameter.  For  example,  the  oscil¬ 
lator  is  centered  at  <w0  =  300  cm-1  (33  cm_1  =  1  THz); 
however,  although  R((o)  can  be  seen  to  start  around  oj0,  it 
asymmetrically  extends  to  frequencies  ~350  cm-1  on  the 
high-6)  end.  As  both  Rico)  and  T(o>)  are  absolute  quanti¬ 
ties,  they  are  constrained  to  lie  between  0  and  1.  Thus  the 
apparent  band  width  is  given  by  a  combination  of  the 
natural  linewidth  (y),  the  oscillator  strength  (>S),  and 
other  contributions  to  the  spectral  weight  at  higher  fre¬ 
quencies  (€,,).  Thus  the  form  of  an  oscillator  in  R(ooi)  and 
T((o)  can  be  quite  complex,  and  the  situation  is  further 
complicated  when  two  oscillators  occur  in  the  same  fre¬ 
quency  range. 

The  reflectance  plotted  in  Fig.  1  is  the  result  of  a 
frequency-dependent  dielectric  contribution,  as  given  by 
Eq.  (2),  at  normal  incidence.  The  situation  is  significantly 
more  complicated  if  frequency-dependent  fi(co)  also  exists 
in  a  material  or  metamaterial.  In  Eqs.  (1)  we  present  the 
Fresnel  equations  that  determine  the  reflectance  from  a 
material  described  by  the  macroscopic  form  of  Maxwells’ 
equations  at  an  angle  8  from  the  surface  normal,  for  both 
transverse  electric  (TE)  and  transverse  magnetic  (TM)  po¬ 
larizations.  As  can  be  seen,  both  e{co)  and  fi(co)  contribute 
to  -R(o)) 


T?(o>)te  - 


cos  0-  /jl  % k/u  -  sin2  8  2 
cos  8  +  /j.~1\[e/ji.  -  sin2  6 


J?(w)tm  - 


e  cos  8  -  yje/j,  -  sin2  8  2 
ecos  (9+  \je/i-  sin2  8 


(1) 


Frequency  (THz)  Frequency  (THz) 


0  150  300  450  0  150  300  450  600 


Frequency  (cm1)  Frequency  (cm1) 

Fig.  1.  Some  optical  constants  for  a  Lorentz  oscillator.  The  os¬ 
cillator  is  modeled  using  Eq.  (2)  with  the  parameters  given  as 
S'  =  50, 000  cm-1,  6oo=  1.5 ,  y- 10  cm-1,  o>0  =  300  cm-1. 
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Frequency  (THz) 


Frequency  (THz) 


Frequency  (cm1)  Frequency  (cnr1) 

Fig.  2.  Oblique  angle  reflectance  plotted  for  0=30°  as  deter¬ 
mined  by  Eqs.  (1)  for  both  TE  and  TM  polarizations.  In  the  left 
panel  the  reflectance  is  a  result  of  a  magnetic  contribution  with 
e=  1  in  this  case.  The  right  panel  displays  B(a>)  for  an  electric  os¬ 
cillator  e(a>)  and  //,=  !. 


At  normal  incidence  (d=0),  the  reflectance  (or  trans¬ 
mission)  spectra  do  not  allow  one  to  discriminate  between 
e(to)  or  /i(oj)  contributions.  That  is,  assuming  magnetic 
and  electric  oscillators  are  of  equal  strength,  the  reflec¬ 
tance  spectra  produced  by  these  two  oscillators  are  iden¬ 
tical.  This  is  why  it  is  important  when  considering 
metamaterials  to  perform  measurements  of  both  the  am¬ 
plitude  and  the  phase  (vector  measurements).  One  pos¬ 
sible  way  to  determine  whether  contributions  are  mag¬ 
netic  or  electric  is  to  perform  ellipsometry.  In  Fig.  2  we 
show  simulated  oblique  angle  reflectance  for  a  magnetic 
oscillator  (left  panel)  and  for  a  electric  resonator  (right 
panel)  both  at  0=  30°.  An  electric  resonance  is  weaker  for 
light  polarization  in  the  TE  mode  compared  with  that  of 
the  TM  mode  and  the  opposite  is  true  for  a  magnetic  reso¬ 
nance.  In  general,  for  an  oblique  incidence  reflectance 
measurement,  the  TE/TM  ratio  for  resonances  in  e(co)  is 
less  than  1  and  greater  than  1  for  modes  associated  with 
/i(o>).  Ellipsometry  inherently  measures  spectra  in  ratios, 
and  thus  it  is  significantly  more  sensitive  than  other 
frequency-domain  spectroscopies. 


3.  ELECTRODYNAMICS  OF 
METAMATERIALS 

A  composite  EM  metamaterial  is  typically  constructed 
from  a  periodic  array  of  one  or  more  conducting  elements 
that  determine  the  overall  EM  response  of  the  entire 
structure.  Straight  wire  or  rod  segments  (left  panel  of  Fig. 
3)  are  used  to  achieve  a  negative  e(co)  response,  while  el¬ 
ements  termed  split-ring  resonators  (SRRs),  right  panel 
of  Fig.  3,  are  often  used  to  obtain  a  negative  /i(oj)  re¬ 
sponse.  Metamaterials  are  microstructured  or  nanostruc- 
tured  materials  designed  to  give  an  effective  EM  response 
that  can  be  described  by  the  macroscopic  form  of  Max¬ 
wells’  equations.  They  are  designed  and  utilized  to  work 
in  the  effective-medium  regime.21  It  is  thus  required  that 
the  individual  elements  that  form  these  composites  be 
much  smaller  than  the  wavelength.  That  is,  if  [a]  is  the 
characteristic  length  of  a  element  from  which  a  metama¬ 
terial  is  formed,  and  A.  is  the  wavelength  at  which  one 
would  like  to  utilize  the  metamaterial,  then  it  is  neces¬ 
sary  that  (a<SA).  This  condition  is  relatively  easy  to  sat¬ 
isfy  at  microwave  frequencies,  and  values  of  A/a —  100 
have  been  demonstrated.  However,  at  higher  frequencies 
the  a<  A  requirement  is  more  difficult  to  achieve  since 


frequency  and  wavelength  scale  inversely.  Typical  values 
at  THz  frequencies  of  A/a  — 10  are  achievable  with  con¬ 
ventional  lithographic  methods. 

An  array  of  wires  or  metal  meshes  is  the  simplest  pos¬ 
sible  metamaterial  that  allows  one  to  obtain  the  regime  of 
an  artificial  dielectric  mentioned  above.22,23  The  work  of 
Ulrich  in  the  late  1960s  and  others  pursued  wire  grid 
meshes  and  their  complementary  structures  to  form  fil¬ 
ters  for  the  IR  or  THz  frequencies.24,25  Other  work  fo¬ 
cused  on  artificial  low-frequency  plasmon  materials,6,26 
which  utilize  thin  straight  conducting  wire  materials  to 
achieve  a  tuned  <op  and  e(co)  response.  Most  efforts  push¬ 
ing  EM  metamaterials  to  higher  frequencies  have  been  fo¬ 
cused  on  the  magnetic  component;  thus  we  will  address 
this  issue  below,  and  only  a  brief  review  of  the  electric 
component  will  be  covered. 

A.  Electric  Metamaterials 

Artificial  wire  composites  are  used  at  frequencies  below 
(op  and  can  be  thought  of  as  a  dilute  plasma.  Essentially 
the  wires,  compared  to  a  bulk  metal,  have  a  reduced  num¬ 
ber  density  in),  as  the  carriers  (electrons)  are  restricted  to 
the  physical  dimensions  of  the  wires.  The  plasma  fre¬ 
quency  is  given  by  o>2  =  4tt -(ne2/m*),  and  the  frequency  de¬ 
pendence  of  the  electric  response  is  given  by  the  Drude- 
Lorentz  model  as 


wp 

e(ai)  =  +  -g - j — : —  >  (2) 

cog  -  (o  -i  y<o 

where  is  the  real  part  of  a  dielectric  constant  at  o> 
=  oo ,  (o0  is  the  center  frequency  of  the  oscillator,  and  y  is 
the  damping.  For  bulk  metals  the  number  density  is  typi¬ 
cally  — 1022  cm3,  and  thus  ojp  usually  occurs  in  the  ultra¬ 
violet.  In  the  case  of  infinitely  long  straight  wire  seg¬ 
ments,  &>o  =  0  and  the  response  is  given  by  the  Drude 
model  (a  Lorentz  oscillator  at  zero  to).27  The  plasma  fre¬ 
quency  of  artificial  wires  in  terms  of  the  geometry  of  the 
media  was  worked  out  in  Refs.  5,23  to  be  <op 
=  (2'7TC2)/(a2  In  a/r),  where  a  is  the  length  of  the  unit  cell, 
r  is  the  radius  of  the  wires,  and  c  is  the  speed  of  light. 
Thus,  for  example,  with  dimensions  such  as  a  =  120  /xm 
and  r-44  /im,  a  plasma  frequency  of  —  1  THz  is  predicted. 

In  the  inset  to  Fig.  4  we  show  an  artificial  straight  wire 
medium  fabricated  for  plasma  frequency  response  in  the 
THz  regime.28  The  radius  of  the  wires  is  15  /xm  and  the 


Fig.  3.  Schematics  of  two  artificial  structures  commonly  used  as 
EM  metamaterials.  Left  panel  shows  a  wire  medium  supporting 
e<  0  for  the  polarization  of  the  incident  EM  wave  as  shown  and 
for  frequencies  io<a>p.  Right  panel  shows  the  unit  cell  for  an  ar¬ 
tificial  magnetic  metamaterial.  The  E  field  lies  along  x ,  B  along 
z,  and  k  along  -y.  By  formation  of  a  bipartite  lattice  utilizing 
both  the  electric  and  the  magnetic  metamaterial,  a  NI  material 
may  be  constructed. 
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Fig.  4.  Transmission  spectra  for  an  electric  metamaterial.  The 
blue  curves  are  for  polarization  perpendicular  to  the  wires,  and 
the  red  is  for  polarization  along  the  wires.  The  transmission 
spectra  near  the  predicted  plasma  frequency  approach  zero,  in¬ 
dicating  screening  of  the  artificial  medium.  Inset  shows  a  picture 
of  the  artificial  wire  medium. 
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Fig.  5.  Drude-Lorentz  model  for  an  oscillator  model  using  Eq. 
(2).  In  (a)  and  (b)  the  real  part  (solid  curve)  and  the  imaginary 
portion  (dashed  curve)  are  given,  (a)  A  Drude  oscillator  (a>0  =  0); 
(b)  the  same  oscillator  but  now  centered  at  o)0=50  cm-1,  (c)  and 
(d)  A  comparison  between  two  magnetic  oscillators  given  from 
Ref.  6  (grey  curves)  and  a  Lorentz  oscillator  (black  curves)  is 
shown.  The  real  portions  are  plotted  in  (c)  and  the  imaginary 
parts  in  (d). 


lattice  spacing  is  120 /xm,  predicting  wp  =  0.7THz.  A 
polarization-dependent  transmission  measurement  is 
shown  in  the  main  panel  of  Fig.  4.  For  light  polarized  with 
the  electric  field  parallel  to  the  wires,  at  frequencies 
above  cop  the  carriers  in  the  wires  are  no  longer  able  to 
screen  the  electric  field  of  an  incident  EM  wave  and  the 
medium  becomes  transparent.  However,  when  the  polar¬ 
ization  state  of  the  light  is  rotated  by  90°  the  medium  is 
transparent,  apart  from  some  dips  in  the  spectrum  re¬ 
lated  to  the  imperfections  in  the  wire  medium.  Thus  the 
measurements  are  in  good  agreement  with  theory.  The 


electric  metamaterial  is  reasonably  in  the  regime  of  the 
effective-medium  limit,  as  \/a  ~  3.6. 

In  Figs.  5(a)  and  5(b)  we  show  the  theoretical  response 
of  both  straight  wire  and  cut  wire  segments.  In  a  continu¬ 
ous  wire  medium,  the  electrons  can  freely  propagate 
along  the  length  of  the  wire  and  can  be  described  by  the 
Drude  model.  In  the  case  of  cut  wires  the  oscillator  in  the 
free  case  is  shifted  to  nonzero  frequencies.  The  cut  wire 
model  is  also  well  described  as,  compared  with  continuous 
wires,  now  having  a  depolarization  field,  which  adds  ca¬ 
pacitance  to  the  model  and  thus  creates  a  resonance.  The 
ability  to  utilize  both  continuous  wires  and  cut  wire  seg¬ 
ments  provides  great  versatility  for  tuning  an  artificial 
electric  metamaterial. 

B.  Magnetic  Metamaterials 

An  element,  the  SRR,  has  been  the  main  component  uti¬ 
lized  to  implement  /x <  0  response.  A  single  SRR  is  de¬ 
picted  in  the  right  panel  of  Fig.  3.  A  SRR  is  a  subwave¬ 
length  resonator,  meaning  that  its  physical  dimensions 
are  much  smaller  than  the  wavelength  at  the  resonant 
frequency  o>0.  An  incident  time-varying  magnetic  field  ap¬ 
plied  as  shown  in  the  right  panel  of  Fig.  3  creates  currents 
circulating  around  the  ring.  The  split  gap  in  the  ring  al¬ 
lows  a  capacitance  to  be  built  up  and  thus  the  resonance 
can  be  thought  of  as  an  inductive-capacitive  (LC)  reso¬ 
nance  with  a)0ocVl/LC.  The  inner  counter-wound  ring 
adds  further  capacitance  and  thus  helps  to  meet  the  effec¬ 
tive  media  requirement.  An  array  of  SRRs  can  be  formed 
and  thus  a  magnetic  metamaterial  can  be  described  by  an 
effective  /x  response. 

The  particular  frequency-dependent  form  of  the  mag¬ 
netic  response  for  SRRs  was  worked  out  in  Ref.  6  to  be 

w2 

=  1  +  F— - - — : — ,  (3) 

ci>q  —  o)  —  lyco 


where  F  is  the  geometric  filling  fraction  of  the  unit  cell, 
related  to  the  dimensions  of  the  SRRs.  This  form  for  a 
resonance  is  not  causal  and  thus  in  this  paper  we  use 

2 

ji((D)  =  1  +  -g - — ,  (4) 

(Oq  —  co  —  i  yco 


where  o\riv  is  the  magnetic  plasma  frequency.  The 
frequency-dependent  magnetic  responses  given  by  Eq.  (4) 
and  that  given  by  Ref.  6  are  quite  similar  [see  Figs.  5(c) 
and  5(d)  for  a  direct  comparison]. 

Although  the  SRR  is  a  relatively  simple  geometric  ob¬ 
ject,  it  is  a  rather  complex  EM  structure.  This  complexity 
stems  from  the  fact  that  these  materials  are  necessarily 
described  by  the  most  general  form  of  the  constitutive  re¬ 
lations.  There  may  be  as  many  as  36  complex  quantities 
to  determine,  and  thus  the  standard  S  parameter  (vector 
reflection  and  transmission  measurements)  yields  incom¬ 
plete  EM  information.  At  THz  and  higher  frequencies 
where  vector  measurements  are  not  common,  retrieval  of 
the  optical  constants  is  even  more  difficult,  although 
methods  such  as  ellipsometry  or  THz  TDS  may  prove  use¬ 
ful.  Typically,  artificial  magnetic  metamaterials  are  con- 
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structed  from  conducting  elements,  and  thus  a  full  EM 
characterization  is  a  necessity,  since  an  electric  response 
is  unavoidable. 

SRRs  have  been  shown  to  be  bianisotropic,29  and  thus  a 
brief  review  of  the  constitutive  relations  for  such  materi¬ 
als  is  in  order.  For  a  complete  review,  the  interested 
reader  is  referred  to  the  work  of  Kong.30  The  constitutive 
relations  for  bianisotropic  materials  may  be  written  as30 
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In  Eq.  (5)  the  permittivities  within  the  2X2  matrix  are 
tensors  of  rank  2.  For  completeness  we  summarize  in 
Table  1  the  various  definitions  for  materials  with  differ¬ 
ent  symmetry  properties  of  the  constitutive  relations.  The 
terms  i;  and  £  are  called  the  magneto-optical  permittivi¬ 
ties,  and  they  describe  coupling  of  the  magnetic-to-electric 
response  and  electric-to-magnetic  response,  respectively. 
Then  for  the  SRRs  with  the  coordinate  system  as  specified 
in  the  right  panel  of  Fig.  3,  and  the  expected  bianisotropic 
nature,  the  material  parameters  are 
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and  the  magneto-optical  permittivities  are 
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The  frequency  dependence  of  Eqs.  (6)  and  (7)  was 
worked  out  to  be29 
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where  the  coefficients  a,  b,  c,  and  d  are  related  to  the  ge¬ 
ometry  of  the  SRR  and  materials  that  are  placed  in  the 
unit  cell.  It  was  further  assumed  that  the  medium  was  re¬ 
ciprocal  and  thus  £yz=-£zy.31  As  mentioned,  since  both  an 
electric  resonance  and  a  magnetic  resonance  can  be  de¬ 
scribed  with  a  Lorentzian  oscillator,  and  moreover  enter 
similarly  in  the  Fresnel  Eqs.  (1)  it  is  difficult  to  distin¬ 
guish  between  the  two  in  either  a  reflection  or  transmis¬ 
sion  experiment  at  normal  incidence. 

Equations  (6)  and  (7)  are  for  SRRs  in  one  dimension; 
however,  the  form  of  the  magneto-optical  permittivities  is 


Table  1.  Material  Definitions  with  Respect  to  the 
Symmetries  of  the  Constitutive  Relations" 


Material  Designation 

6,  fj. 

4  C 

Isotropic 

I 

1 1=0 

Anisotropic 

e  and/or  Jjl  not  «:I 

1 1=0 

Biisotropic 

ocl 

cel 

Bianisotropic 

All  other  cases 

All  other  cases 

a\  denotes  the  identity  matrix. 


significantly  more  complicated  in  three  dimensions  (3D). 
Although  the  magnetic  response  generalizes  to  3D  and  is 
isotropic,  the  bianisotropic  SRRs  add  two  terms  to  the 
magneto-optical  permittivities,  e.g.,  4 y=~lyz  and  Cyx=lx y> 
and  another  frequency-dependent  electrical  response, 
e.g.,  ezz(co),  depending  upon  how  the  gap  in  the  SRR  is  ori¬ 
ented.  Thus  the  SRR  is  not  an  ideal  magnetic  material  to 
extend  to  higher  dimensions. 


4.  SPECTROSCOPY  OF  METAMATERIALS 

The  current  effort  to  extend  metamaterials  to  higher  fre¬ 
quencies  requires  innovative  approaches  to  characteriza¬ 
tion.  Characterization  methods  that  require  the  direct 
measurement  of  the  transmitted  and  reflected  phases  are 
less  convenient  at  submillimeter  wavelengths.  To  circum¬ 
vent  this  limitation,  many  researchers  have  made  use 
solely  of  transmittance  or  reflectance  measurements 
coupled  with  numerical  simulations  to  characterize  the 
response  of  high-frequency  metamaterials. 

The  characterization  of  high-frequency  metamaterials 
is  further  complicated  by  the  constraints  of  fabrication.  To 
date,  fabrication  techniques  have  largely  been  based  on 
planar  lithography,  in  which  a  single  layer  of  a  metama¬ 
terial  structure  is  deposited  on  a  substrate.  For  limited 
quantity  proof-of-concept  experiments,  the  processing  is 
expensive  and  time-consuming,  such  that  one  plane  of  the 
sample  is  often  all  that  is  feasible  to  produce.  The  assign¬ 
ment  of  bulk  material  parameters  to  a  single  plane  of  a 
sample  adds  a  substantial  burden  to  the  characterization 
effort. 

In  the  absence  of  phase  information,  ellipsometry  is  a 
promising  approach  to  material  characterization  at  high 
frequencies.  Yen  et  al.  have  performed  a  partial  ellipso- 
metric  study  of  a  planar  SRR  sample  designed  to  operate 
at  THz  frequencies.  This  experiment  illustrates  both  the 
possibilities  as  well  as  the  difficulties  associated  with  el- 
lipsometric  analysis  of  complex  materials.11  Three  differ¬ 
ent  samples  of  3  /jcm  thick  copper  SRRs  were  fabricated 
on  a  400  /r m  thick  quartz  substrate  (inset  to  Fig.  6)  and 
each  were  characterized  with  both  TE-  and  TM-polarized 
light.  Additionally,  since  SRRs  are  known  to  be  bianisotro¬ 
pic,  the  samples  were  rotated  by  90  deg  in  the  plane  and 
were  again  characterized  with  both  polarizations.  The  re¬ 
sults  of  measurements  of  structures  with  different  dimen¬ 
sions  and  geometries  directly  demonstrate  the  scalability 
of  SRRs  over  a  band  of  THz  frequencies,  as  shown  in 
Fig.  6. 
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Fig.  6.  Measurements  of  micrometer-sized  SRRs  designed  for 
magnetic  response  at  THz  frequencies  after  Ref.  11.  Inset  to  bot¬ 
tom  panel  shows  a  picture  of  one  SRR,  and  the  corresponding  di¬ 
mensions:  gap  between  the  inner  and  outer  ring  (G),  the  width  of 
the  metal  lines  (W),  the  length  of  the  outer  ring  (L),  and  the  lat¬ 
tice  parameters  were  2,  4,  26,  and  36  /im  for  the  sample  exhibit¬ 
ing  the  highest  frequency  resonance  (red  curve);  3,  4,  32,  and 
44  fj,m  (black  curve),  and  3,  6,  36,  and  50  /im  for  the  sample  with 
the  lowest  resonant  frequency  (blue  curve).  Top  panel  displays 
the  ellipsometric  parameter  described  in  the  text  for  the  three 
SRR  samples  (red,  black,  and  blue  curves).  The  dashed  black 
curve  and  the  solid  gray  curve  are  for  the  same  measured  param¬ 
eter  for  copper  and  quartz,  respectively.  Middle  and  bottom  pan¬ 
els  display  the  real  and  imaginary  portions  of  /z(a>),  respectively, 
as  simulated  for  the  three  different  SRRs. 

The  thin  SRR  planar  geometry  by  definition  is  aniso¬ 
tropic  and  of  finite  thickness,  so  that  the  Fresnel  scatter¬ 
ing  formulas  become  considerably  more  complicated  than 
those  in  Eqs.  (1).  Moreover,  because  the  planar  sample 
does  not  possess  a  plane  of  symmetry,  the  magnetoelectric 
coupling  term  cannot  be  neglected.  With  at  least  six  un¬ 
known  significant  material  parameters  (real  and  com¬ 
plex)  to  be  determined,  transmittance  or  reflectance  mea¬ 
surements  on  the  structure  at  several  angles  of  incidence 
would  be  required  to  properly  characterize  the  structure. 

The  ellipsometry  ratio  (the  ratio  of  TM-to-TE  polarized 
light),  plotted  in  the  top  panel  of  Fig.  6,  demonstrates  that 
the  SRRs  have  an  enhanced  response  at  the  predicted 
THz  frequencies,  consistent  with  simulations  on  the  iden¬ 
tical  geometries.  By  correlating  the  observed  resonant  fre¬ 
quency  with  that  predicted  from  simulations,  a  strong 
case  for  the  magnetic  response  of  SRRs  can  be  made. 
Given  the  complicated  manner  in  which  all  the  resonant 
terms  enter  into  the  ellipsometry  formulas,  a  more  quan¬ 
titative  statement  about  the  permeability  component  can¬ 
not  be  made.  However,  the  magneto-optical  terms  (£  and 
Q  manifest  themselves  in  spectroscopy  as  polarization  ro¬ 
tation  or  mixing,  resulting  in  elliptically  polarized  light. 


Signs  of  this  cross  polarization  were  not  observable 
through  cross-polarized  measurements. 

Although  the  bianostropic  planar  samples  studied  in 
Ref.  11  lead  to  a  complication  in  characterization,  the 
situation  can  be  greatly  improved  by  fabricating  samples 
that  possess  symmetry  planes  in  all  directions.  By  analyz¬ 
ing  the  reflectance  of  symmetric  structures  for  waves  in¬ 
cident  at  nonnormal  incidence,  and  with  careful  consider¬ 
ation  of  the  incident  polarization  of  the  wave  relative  to 
the  structure,  a  quantitative  retrieval  of  the  permeability 
should  be  quite  feasible. 

By  applying  ever-improving  microfabrication  and  nano¬ 
fabrication  methods,  researchers  have  made  steady 
progress  in  advancing  the  frequency  response  of  metama¬ 
terials  upwards  in  frequency  over  the  past  several  years. 
A  study  in  2004  (Ref.  32)  reported  a  far-infrared  reso¬ 
nance  with  aluminum  SRRs  0.5  /xm  thick  on  a  0.7  /im 
thick  silicon  substrate  at  ~28  THz.  The  characterization 
of  the  SRRs  was  based  on  a  transmission  experiment  with 
several  laser  lines  X.j  =  10.86  |iim,  \2  =  9.43  /im,  and  \3 
=  808  nm.  The  characteristic  length  of  the  SRRs  was  a 
=  10.48  fim  and  the  lattice  parameter  was  ~15  /mi .  The 
measurements  from  Ref.  32  indicate  that  there  is  a  differ¬ 
ence  in  absorption  between  the  SRR  orientation  for  which 
a  magnetic  resonance  is  expected  compared  with  the  SRR 
orientation  for  which  no  resonance  is  expected.  Although 
the  SRR  samples  used  in  Ref.  32  are  strictly  outside  of  the 
effective-medium  regime  (since  wavelengths  \1;  \2,  and  X3 
yield  a  ratio  of  Kl!a ~  1),  the  scattering  results  are  consis¬ 
tent  with  an  enhanced  resonant  response  in  the  SRR  over 
the  20-30  THz  frequency  range.  By  reducing  the  size  of 
the  unit  cell  relative  to  the  wavelength  and  symmetrizing 
the  SRR  structures,  transmittance  measurements  com¬ 
bined  with  various  sample  rotations  could  result  in  a 
novel  characterization  method. 

Shortly  after  the  measurements  of  Chuan  et  al.32  were 
presented,  Linden  et  al.  made  use  of  the  inherent  magne¬ 
toelectric  coupling  in  SRRs  to  characterize  a  SRR 
metamaterial  at  100  THz.12  In  this  work,  transmittance 
and  reflectance  measurements  were  carried  out  at  near¬ 
normal  incidence  for  three  different  planar  SRR  samples, 
all  having  the  same  dimensions  (a~320nm),  but  with 
varying  lattice  parameters.  In  these  experiments  the 
magnetic  field  was  oriented  so  as  to  lie  in  the  plane  of  the 
SRR  arrays  and  the  electric  field  was  oriented  along  ei¬ 
ther  the  asymmetric  direction  y  axis  [left  panel  of  Fig. 
7(c)]  or  the  symmetric  x  axis  [right  panel  of  Fig.  7(c)].  Be¬ 
cause  the  incident  magnetic  field  is  perpendicular  to  the 
axes  of  the  SRRs,  the  magnetic  resonance  of  the  SRR  is 
not  driven  directly,  and  the  magnetic  response  of  the  SRR 
must  be  assessed  indirectly  by  magnetoelectric  coupling 
as  follows.  When  the  polarization  of  the  electric  field  is  di¬ 
rected  along  the  asymmetric  y  axis,  a  peak  in  R(ooi)  is  re¬ 
corded.  However,  when  the  polarization  is  directed  along 
the  symmetric  direction,  this  peak  disappears.  The  peak 
in  R(co)  and  the  corresponding  dip  in  T(a>)  clearly  reveal 
the  expected  SRR  resonance  that  occurs  at  a  wavelength 
near  3  /zm.  As  a  control,  the  authors  also  measure  closed 
rings  of  otherwise  identical  dimensions  and  demonstrate 
the  absence  of  a  resonant  response  over  the  same  wave¬ 
lengths.  The  resonance  of  the  SRR  samples  is  consistent 
with  simulations  performed  by  the  authors,  which  reveal 
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that  rotating  currents — a  prerequisite  for  artificial 
magnetism — are  excited  within  the  SRRs. 

The  transmittance  experiments  of  Ref.  12  provide  di¬ 
rect  evidence  of  enhanced  SRR  response  at  wavelengths 
nearing  the  visible.  Yet,  because  of  the  complicated  nature 
of  the  bianosotropic  SRR  sample,  which  possesses  both 
resonant  permittivity  as  well  as  resonant  permeability 
and  magnetoelectric  coupling  [see  Eqs.  (7)],  it  is  difficult 
to  extricate  quantitatively  the  desired  magnetic  response. 
Nevertheless,  the  correlation  of  the  SRR  resonance  with 
rotating  currents  implies  that  artificial  magnetic  struc¬ 
tures  can  be  designed  and  fabricated  to  operate  at  wave¬ 
lengths  much  smaller  than  were  first  deemed  possible.  Al¬ 
though  the  expected  magnetic  response  was  shown  to  be 
considerably  weaker  than  the  electric  response  in  these 
particular  samples,  altering  the  symmetry  of  the  SRR  de¬ 
sign  could  result  in  a  material  whose  electric  and  magne¬ 
toelectric  responses  are  suppressed. 

Utilizing  a  novel  fabrication  process  based  on  LIGA, 
Moser  et  al.  formed  a  metamaterial  composed  of  planar 
SRRs  and  discontinuous  wires.33  The  authors  observed  a 
resonant  electromagnetic  response  at  ~2.5  THz.  Several 
different  metamaterials  were  fabricated  using  both  Ni 
and  Au  to  form  SRR  or  wire  elements.  The  typical  dimen¬ 
sions  of  the  SRRs  were  ~80  jj.m.  All  the  structures  were 
embedded  within  a  photoresist,  transparent  over  the  THz 
band,  and  characterized  with  Fourier-transform  IR  trans¬ 
mission  measurements.  Simulations  were  performed  on 
the  composite  structures  and  compared  with  analytical 
predictions  using  the  formulas  of  Pendry  et  al.5  The  ex¬ 
perimental  transmission  measurements  of  the  composite 
metamaterials  revealed  a  peak  in  T(io),  in  good  agree- 
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Fig.  7.  SRR  arrays  designed  to  operate  in  MIR  and  near-IR  fre¬ 
quencies  (after  Ref.  12).  (a)  A  photograph  of  a  typical  SRR  cell 
and  (b)  of  an  array  used  for  spectroscopic  studies  in  MIR  near-IR 
frequencies,  (c)  The  transmission  (red)  and  reflectance  (blue) 
spectra.  The  two  polarization  configurations  are  shown  on  top  of 
the  two  columns.  Lattice  constant  a  =  450nm  corresponds  to 
nominally  identical  SRRs. 
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Fig.  8.  Magnetic  metamaterial  constructed  from  Au  designed 
for  resonance  at  MIR  frequencies.13  (a)  Depicts  how  charges  will 
flow  in  the  device  upon  the  application  of  a  time-varying  mag¬ 
netic  field  directed  along  the  z  axis,  (b)  Shows  the  currents  from 
(a)  as  well  as  the  resulting  image  charges.  The  charges  and  im¬ 
age  charges  thus  form  a  net  circulating  current  resulting  in  a 
magnetic  response.  Upon  excitation  by  an  external  electric  field 
E  directed  along  they  axis,  charges  will  flow  as  shown  (c),  as  well 
as  image  charges  shown  in  (d),  and  thus  a  net  circulation  and 
magnetic  field  B  will  be  created. 


ment  with  that  predicted  by  simulation  and  the  analytic 
formulas.  The  structures  were  simulated  with  polarized 
light,  but  the  experiment  was  carried  out  with  unpolar¬ 
ized  light.  Characterization  of  the  composite  metamateri¬ 
als  was  completed  with  the  light  propagating  along  the 
surface  normal  of  the  structures,  and  thus  the  magnetic 
field  was  in  plane.  As  the  experiment  was  performed  with 
unpolarized  light,  and  due  to  the  orientation  of  the 
sample  with  respect  to  the  incoming  beam,  contributions 
to  T(a j)  were  due  to  both  the  eyy  and  £zy  terms  of  Eq.  (5). 

The  examples  above  illustrate  both  the  successes  as 
well  as  the  difficulties  in  scaling  artificial  magnetic  and 
NIMs  to  high  frequencies.  As  we  have  seen,  one  of  the 
more  complicating  factors  is  the  inherent  bianisotropy  of 
the  asymmetric  SRR  unit  cell,  which  places  a  greater  bur¬ 
den  on  retrieval  techniques.  Away  to  avoid  the  complica¬ 
tion  of  bianisotropy  is  to  symmetrize  either  the  individual 
SRRs  or  arrange  asymmetric  SRRs  such  that  they  are 
symmetric  on  average.34  SRRs  designed  to  minimize  mag¬ 
netoelectric  coupling  at  IR  frequencies  have  been  re¬ 
ported,  but  have  yet  to  be  measured.35  Structures  that 
bear  some  similarity  to  those  suggested  in  Ref.  35  have 
been  fabricated  and  measured  at  MIR  frequencies.13  The 
staplelike  element  that  forms  the  basis  of  these  IR  mate- 
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rials  is  depicted  in  Fig.  8.  The  staple  utilizes  a  ground 
plane  to  complete  the  other  half  of  the  structure,  thus 
achieving  an  effective  magnetic  response.  These  struc¬ 
tures  have  a  characteristic  length  of  ~200  nm  and  a  pre¬ 
dicted  resonant  frequency  of  ~60  THz.  The  Au  staples  lie 
upon  a  dielectric  material  of  80  nm  thickness  above  a 
ground  plane.  The  structure  is  excited  by  an  incident  EM 
wave  with  the  magnetic  component  lying  along  the  z  axis 
as  shown  in  Fig.  8(a).  The  time-changing  magnetic  field 
causes  circulation  of  charges  in  the  structure.  The  ground 
plane  allows  charges  (currents)  to  run  in  the  image  plane 
in  the  opposite  direction,  thus  resulting  in  a  net  circula¬ 
tion  of  charges  and  yielding  a  magnetic  response. 

In  Figs.  9(a)  and  9(b)  we  show  pictures  of  some  of  the 
magnetic  structures  of  Ref.  13.  Figure  9(c)  shows  reflec¬ 
tion  spectra  for  one  sample  (solid  curve)  and  simulated  re¬ 
flectance  (dotted  curve).  The  response  of  one-dimensional 
structures  was  simulated  using  rigorous  coupled-wave 
analysis  (RCWA),  and  bulk  properties  of  Au  were  used. 
The  extracted  frequency-dependent  /x(<u)  is  shown  in  Fig. 
9(d).  Good  agreement  between  the  expected  resonance 
and  the  experiment  is  apparent. 

Interestingly,  although  upon  first  inspection  the  staple 
structure  would  appear  to  be  a  symmetric  structure  and 
hence  free  of  bianisotropy,  because  of  the  asymmetric 
drive  introduced  by  the  ground  plane,  the  staple  structure 
actually  exhibits  the  same  magneto-optical  response  as 
the  asymmetric  SRRs.  In  Fig.  8(c)  it  is  shown  how  charges 
will  flow  under  the  application  of  an  electric  field.  Image 
charges  will  flow  in  the  opposite  direction,  and  a  net  cir¬ 
culation  of  charges  will  result  and  a  magnetic  response 
will  be  exhibited.  Thus  just  as  the  SRRs  are  bianisotropic, 
so  is  the  staple  structure.  This  turns  out  to  be  of  little 
practical  consequence  in  this  particular  experiment,  how¬ 
ever,  since  the  driving  electric  field  is  substantially  re¬ 


duced  by  the  presence  of  the  ground  plane.  Thus,  al¬ 
though  the  staple  structure  exhibits  bianisotropy,  the 
reduction  of  the  driving  electric  field  enables  the  magnetic 
properties  of  the  SRR  to  be  probed  more  or  less  directly. 
Simulations  on  the  staple  structure  suggest  that  a  fre¬ 
quency  region  in  the  IR  should  exist  over  which  the  effec¬ 
tive  permeability  is  negative. 

The  above-discussed  measurements  of  metamaterials 
have  involved  variations  of  the  experimental  techniques 
(i)  and  (iii)  discussed  in  Section  2.  As  mentioned,  other 
techniques  such  as  interferometric  methods  (iv)  or  THz 
TDS  (v),  while  less  available  and  convenient,  may  provide 
a  more  complete  material  characterization.  These  meth¬ 
ods  benefit  by  the  addition  of  phase  information.  Recently, 
measurements  of  materials  designed  to  exhibit  NI  at  X 
~2.0  fini  were  presented  using  interferometry  in  a 
near-IR  range.36  These  metamaterials  are  different  from 
typical  metamaterials  used  at  microwave  and  THz  fre¬ 
quencies,  consisting  of  a  square  array  of  through  holes  in 
two  conducting  layers  separated  by  a  dielectric  all  on  a 
glass  substrate  [see  Fig.  10(a)].  The  lattice  parameter  is 
a  =  838  nm  and  the  diameter  of  the  holes  is  c/  =  360  nm. 
The  conducting  Au  layers  are  30  nm  thick  and  the  dielec¬ 
tric  is  made  of  A1203  and  is  60  nm  thick.  Interferometric 
information  was  obtained  by  forming  two  additional 
samples,  one  by  forming  wires  of  the  metamaterial  on  the 
glass  substrate  and  the  other  by  adding  Au  wires  with  a 
particular  spacing  on  top  of  the  metamaterial.  Then  both 
T(co)  and  R(co)  were  characterized  for  all  samples.  Apeak 
in  T(to)  for  the  metamaterial  sample  occurring  near  2  fim 
is  changed  to  be  a  dip  for  the  metamaterial  wire  sample. 
This  is  taken  to  indicate  a  large  difference  in  phase  be¬ 
tween  the  two  samples.  Both  T(co)  and  R(co)  were  simu¬ 
lated  for  the  samples  using  RCWA,  and  good  agreement 
with  the  experiment  was  obtained.  Finally,  extraction  of 
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Fig.  9.  Magnetic  metamaterial  designed  for  response  at  MIR  frequencies.  Two  different  structures  are  shown  in  (a)  and  (b),  all  con¬ 
structed  from  Au.  In  (c)  reflection  measurements  are  shown  for  sample  c,  not  shown.  In  (d)  the  calculated  real  and  imaginary  parts  of 
fi(a>)  are  shown. 
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(a) 


(b) 


Fig.  10.  Metamaterials  used  for  response  at  near-IR  and  optical 
frequencies,  (a)  A  schematic  of  the  mulitlayer  structure  used  in 
Ref.  36  consisting  of  an  A1203  dielectric  layer  between  two  Au 
films  with  a  square  array  of  holes  (pitch  of  838  nm,  diameter  of 
360  nm).  (b)  Shows  the  pairs  of  Au  nanorods  used  in  Ref.  19  sepa¬ 
rated  by  a  Si02  layer  on  top  of  a  glass  substrate  separated  by  a 
180  nm  layer  of  indium  tin  oxide  (ITO). 

the  index  ( n )  and  impedance  (z)  was  obtained  from  the  ex¬ 
perimental  data  and  it  was  found  to  exhibit  NI. 

Another  phase-dependent  study  focused  on  pairs  of  Au 
nanorods  separated  by  a  dielectric  material,  engineered 
for  NI  response  near  1.2  /mi. 19  The  length  and  width  of 
the  nanorods  was  750  and  170  nm,  respectively,  and  had  a 
thickness  of  50  nm  [see  Fig.  10(b)].  A  metamaterial  was 
formed  by  fabricating  two  nanorods  separated  by  50  nm  of 
a  Si02  on  180  nm  of  indium  tin  oxide  on  a  glass  substrate. 
T(to)  and  R{to)  were  measured  with  a  grating  spectrom¬ 
eter  using  polarized  light,  while  the  phase  was  obtained 
by  using  polarization  and  walk-off  interferometers.  The 
samples  were  simulated  using  finite-difference  time- 
domain  methods,  and  good  agreement  with  the  experi¬ 
ment  was  obtained.  Inversion  of  complex  T(io)  and  R(coi) 
data  indicated  that  a  NI  was  obtained  with  a  maximum 
value  of  zi  =  -0.3  at  X  ~  1.5  /cm. 

5.  TUNABLE  NEGATIVE  INDEX  AND 
METAMATERIALS 

So  far  the  practical  implementation  of  NI  phenomena  has 
been  successfully  accomplished  only  in  microstructured 
and  nanostructured  systems  described  above.  The  re¬ 
sponse  of  a  given  metamaterial  structure  is  therefore  rig¬ 
idly  determined  by  the  geometric  dimensions  of  the  split 
rings  and  other  elements.  From  the  viewpoint  of  applica¬ 
tions  as  well  as  for  the  purpose  of  the  in-depth  analysis  of 
the  electrodynamics  of  metamaterials,  it  is  desirable  to 
implement  structures  with  tunable  resonant  properties. 
In  this  section  we  will  briefly  outline  some  of  the  possible 
solutions  for  the  realization  of  tunable  NIM  systems. 

We  start  by  exploring  the  expression  for  the  resonant 
frequency  of  a  SRR  array6: 


VttVo Cr3’  (9) 

where  r  is  the  ring  radius,  /t0  is  the  static  magnetic  per¬ 
meability,  and  C  is  the  capacitance  per  unit  area  between 
the  two  rings.  While  this  analytic  expression  is  only  ap¬ 
proximate,  detailed  simulations  have  confirmed  that  the 
resonant  frequency  of  the  SRR  is  indeed  sensitive  to  this 
capacitance,  which  in  turn  depends  on  the  value  of  the  di¬ 
electric  constant  of  the  substrate  es.3t  Therefore  the  net 
result  of  a  substrate  dielectric  that  deviates  form  its  free- 
space  value  is  to  renormalize  the  frequency  position  of  the 
resonance  by  the  factor  °cl /\[es.  Therefore  the  easiest 
route  to  achieve  tunable  resonant  behavior  is  through  a 
modification  of  the  dielectric  constant  of  a  substrate.  SRR 
arrays  and  other  types  of  metamaterials  can  be  incorpo¬ 
rated  in  metal-insulator-semiconductor  (MIS)  structures. 
Applying  a  dc  electric  field  between  the  ring  arrays  and  a 
semiconducting  substrate  allows  tuning  the  dielectric 
constant  of  the  insulating  layer  in  a  MIS  device  provided 
that  the  insulator  is  fabricated  from  a  high  dielectric  con¬ 
stant  or  ferroelectric  material.38  To  maintain  the  electric 
field  across  the  insulator  in  the  area  within  the  splits,  it  is 
desirable  to  fill  these  structures  with  semiconducting 
polymers.  Charge  injection  into  a  polymer  allows  one  to 
achieve  an  electric  field  in  MIS  structures  over  large  ar¬ 
eas  without  appreciably  changing  the  dc  conductivity  of 
the  polymer.38  Therefore  in  MIS-based  devices  the  dielec¬ 
tric  constant  and  therefore  the  resonant  frequency  can  be 
tuned  by  varying  the  applied  dc  voltage.  One  can  easily 
envision  that  each  of  the  rings  in  a  SRR  array  may  be  bi¬ 
ased  with  a  different  voltage.  This  selective  tuning  of  the 
dielectric  constant  may  allow  one  to  achieve  the  con¬ 
trolled  distribution  of  the  resonant  frequencies  over  a  pla¬ 
nar  array. 

Alternatively,  the  dielectric  constant  of  the  substrate 
can  be  modified  by  photodoping,  thus  enabling  a  dynami¬ 
cal  modification  of  metamaterial  properties  in  accord  with 
Eq.  (9).  Mode-locked  lasers  readily  allow  one  to  achieve 
high  densities  of  photogenerated  carriers  in  common 
wide-gap  insulators,  over  subpicosecond  time  scales.39  In 
fact,  in  Ref.  39  the  authors  succeeded  in  generating  me¬ 
tallic  Drude  absorption  that  should  be  sufficient  to  dy¬ 
namically  shunt  the  gap  in  the  SRRs.  Therefore  short  la¬ 
ser  pulses  can  be  exploited  for  switching  off  the  magnetic 
response  of  a  SRR  ring.  THz  TDS  is  ideally  suited  to 
probe  the  optical  constants  of  metamaterials  under  pho¬ 
toexcitation.  The  structures  of  Ref.  13  pictured  in  Figs.  8 
and  9  allow  for  the  easy  incorporation  of  dielectrics  and 
other  natural  materials  within  the  gaps  of  magnetic 
metamaterials.  The  added  benefit  of  the  ground  permits 
dc  bias  of  materials  between  the  gaps. 


6.  SUMMARY  AND  OUTLOOK 

We  have  detailed  the  expected  spectroscopic  response 
from  metamaterials  consisting  of  both  electric  and  mag¬ 
netic  artificial  materials.  The  frequency-dependent  form 
for  each  of  these  structures  has  been  analyzed  and  the  dif¬ 
ficulty  of  retrieval  of  the  optical  constants  from  experi¬ 
mental  observables  has  been  stressed.  It  is  suggested  that 
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future  characterization  of  metamaterials  at  THz  and 
higher  frequencies  take  advantage  of  vector  measure¬ 
ments  (phase  and  amplitude),  such  as  ellipsometry,  THz 
TDS,  or  Mach-Zehnder  interferometry.  The  artificial 
magnetic  metamaterial  in  particular  is  a  difficult  struc¬ 
ture  to  thoroughly  characterize.  This  is  due  to  the  fact 
that,  although  it  is  designed  to  give  a  magnetic  resonant 
response  at  a  frequency  &>0,  it  exhibits  an  undesirable 
electric  resonance  at  the  same  frequency,  thus  complicat¬ 
ing  its  use  and  characterization.  A  simpler  alternative 
route  is  use  of  more  symmetric  structures  that  do  not  ex¬ 
hibit  a  bianisotropic  response. 

The  extension  of  metamaterials  to  the  dynamical  re¬ 
gime  has  been  reviewed.  Several  possible  methods  have 
been  suggested,  all  of  which  are  particularly  suited  for 
THz  frequencies,  a  region  of  the  EM  spectrum  lacking  sig¬ 
nificant  natural  material  response. 
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